Small-Scale Cosmic Microwave Background Polarization from Reionization 
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We discuss fluctuations in the cosmic microwave background (CMB) polarization due to scattering 
from reionized gas at low redshifts. Polarization is produced by re-scattering of the primordial 
temperature anisotropy quadrupole and of the kinematic quadrupole that arises from gas motion 
transverse to the line of sight. We show that both effects produce equal E- and B-parity polarization, 
and are, in general, several orders of magnitude below the dominant polarization contributions 
at the last scattering surface to E-modes or the gravitational-lensing contribution to B-modes at 
intermediate redshifts. These effects are also several orders of magnitude below the B polarization 
due to lensing even after subtraction with higher-order correlations, and are thus too small to 
constitute a background for searches for the polarization signature of inflationary gravitational 
waves. 



I. INTRODUCTION 

The angular power spectrum of cosmic microwave 
background (CMB) temperature fluctuations is now be- 
coming a powerful cosmological probe [Jj], both due to 
our detailed understanding of physics during the recom- 
bination era Q and progress on the experimental front 
PL In addition to the primary anisotropies generated 
at the last-scattering surface, CMB photons are also af- 
fected by large-scale structure at low redshifts. These lat- 
ter contributions result from scattering off free electrons 
in clusters or the reionized IGM and from modifications 
due the evolving gravitational field associated with the 
formation of structures 0] . 

The existence of such secondary signals has now be- 
come evident with the recent detection of small-scale 
power in excess of that from primary fluctuations [^J. 
The simplest and most plausible explanation for this 
small-scale power is the Sunyaev-Zel'dovich (SZ; Q), re- 
scattering of CMB photons from hot electrons in un- 
resolved galaxy clusters. Although the power observed 
is considerably larger than theoretical expectations, the 
excess can be accommodated with a relatively small in- 
crease in the power-spectrum amplitude [0 [10[ . 

In addition to small-scale anisotropies in the tempera- 
ture, increasing attention is now being devoted to detec- 
tion of the CMB polarization. Besides resolving cosmo- 
logical parameter degeneracies 0, the polarization will 
allow several unique cosmological and astrophysical stud- 
ies to be carried out. These include a reionization signal 
PI , probes of gravitational lensing |l2|, and a signa- 
ture of inflationary gravitational waves (IGW) through 
its contribution to the B, or curl, modes of the polariza- 
tion 0. 

Given the rapid pace of experimental progress and the 
ride of thumb that the polarization is typically 10% of the 
temperature anisotropy, it is appropriate to investigate 
the polarization produced by the secondary effects that 
have produced the recently detected small-scale power. 
This polarization is produced by Thomson scattering of 



the quadrupole moment of the radiation incident on the 
scatterer. The quadrupole moment can be either the pri- 
mordial quadrupole that the scatterer sees @ or the 
kinematic quadrupole that arises from quadratic terms 
in the Doppler shift when the gas moves transverse to 
the line of sight ||. Small-scale angular fluctuations in 
this polarization are produced by variations in the optical 
depth as a function of position across the sky. 

Since the secondary polarization signals can affect cos- 
mological studies involving the primordial polarization, 
and, by themselves, may provide important information 
on astrophysics at late times, it is important that we both 
quantify and understand the extent to which large-scale 
structure can be a potential source of CMB polarization. 
Here, we discuss the angular power spectra for scattering 
from reionized electrons and study how these may affect 
potentially interesting studies with CMB polarization. 

Our calculation parallels that of Ref. Jig ], although 
differs in that we present a simplified derivation of the 
results based on a flat-sky approximation and use a halo- 
clustering approach to describe fluctuations in the 
electron density, in analogy to similar recent calculations 
of small-scale temperature fluctuations from unresolved 
clusters |||l7],|l8|]. We also include for the first time the 
frequency dependence in the small-scale polarization us- 
ing results for the polarization from individual clusters 

EMU- 

The paper is organized as follows. In § O, we intro- 
duce polarization signals associated with the scattering 
of the primordial CMB temperature quadrupole and with 
the kinematic quadrupole generated by transverse mo- 
tions. We then display and discuss our results in § III. 
Though we present a general discussion of the polar- 
ization, when illustrating results, we will use the cur- 
rently favored ACDM cosmology with matter density (in 
units of the critical density) Q m — 0.35, baryon density 
Qb = 0.05, vacuum-energy density fl\ = 0.65, Hubble 
constant (in units of 100 km sec -1 Mpc -1 ) h = 0.65, 
and spectral n = 1 for primordial density perturbations. 
We employ natural units with the speed of light c = 1 
throughout. 
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II. POLARIZATION POWER SPECTRA 

If the radiation incident on a rcionized electron has a 
nonzero quadrupole moment, then the scattered radia- 
tion will be linearly polarized. The two dominant origins 
for this quadrupole moment are: (a) A quadrupole from 
primordial CMB fluctuations, and (b) a quadrupole from 
the quadratic term in the Doppler shift if the scattering 
gas has a transverse velocity. 

The statistics of the polarization produced in this way 
is inherently non-Gaussian: Although the amplitude of 
the polarization in a given direction is determined by the 
optical depth to Thomson scattering, the orientation is 
determined by the radiation quadrupole incident on the 
scattering gas. The power spectrum is thus in principle 
determined by a convolution of the optical-depth field 
with the quadrupole field. 

In practice, however, we simplify by assuming that 
the quadrupole moment is smooth over large distances, 
and that there are small-scale variations only in the free- 
electron density. We can then proceed to calculate the 
polarization fluctuations that are linear in the electron- 
density fluctuations. The correlation length of the CMB 
quadrupole is comparable to the horizon, and so our as- 
sumptions are fully justified for this case, as long as we 
restrict our attention to multipole moments / > 10. On 
the other hand, the peculiar velocity is correlated only on 
much smaller scale, ~60 Mpc in our fiducial ACDM cos- 
mology. Our calculations of the polarization power spec- 
trum for the kinematic quadrupole will thus be reliable 
only for scales smaller than this, or multipole moments 
I > 100. 

We first discuss the power spectra. A polarization map 
will consist of a measurement of the Stokes parameters 
Q(h) and U(fi) as functions of position n on some patch 
of the sky. We can construct Fourier components Q(l) 
and £7(1) by 

X(Y) = J d 2 h e ll fl A(n), (1) 

where X = Q,U, and 1 is a vector in the plane of a region 
of the sky sufficiently small to be considered flat. Since 
the Stokes parameters Q and U depend on the choice of 
axes, we consider the rotationally invariant combinations, 

E(l) = cos(2<fr)Q(l) + sin(2<fr)Z7(l) 

B(l) = cos(20i)C/(l) - sin(20 1 )Q(l) , (2) 

where <j)\ is the angle between 1 and the chosen x-axis in 
the plane of the sky. We define the angular power spectra 
Cf E and Cf B from 

{Y{\)Y{V)) = {2-Kf5 D {\ + V)Cr , (3) 

where Y = E, B, and the angle brackets denote an aver- 
age over all realizations of the density field. 



As discussed above, we suppose that the radiation 
quadrupole moment is smooth over the region of sky we 
are considering. If so, then 

(E(l)E(V)) = (B(l)B(V)) 

= i((Q(l)Q(l')) + Ml)t/(l'))) , (4) 

while 

(E(l)B(V)) = \ «Q(l)tf(l')> - <Q( W)» = , (5) 

where the latter equality is consistent with parity con- 
servation. These results can be derived by noting that if 
the quadrupole moment is constant, then the orientation 
of the polarization is constant. If so, we may choose our 
axes on the sky so that [7 = 0. Then, 75(1) oc sin(2</>i), 
and 73(1) oc cos(20i), but when averaged over the orien- 
tation angle <f>\, we recover Eq. (^). 

We now proceed to calculate the power spectra induced 
by reionization. The polarization in direction n due to 
scattering from free electrons is an integral along the line 
of sight 

Q(n) - iU(n) = ^-J dr^^a 22 (r) , (6) 

where r is the comoving distance, (dr / dr)(rxi,r) = 
<JTn e (rii, r)a(r), a(r) is the scale factor at a comoving 
distance r, n e (rn,r) is the free-electron density at di- 
rection n at distance r, and ot is the Thomson cross 
section. 

In Eq. (^) , (r) is the radiation quadrupole moment 
at distance r. More precisely, 022 (r) is the coefficient of 
the spherical harmonic Y 22 ifi, 4>) m a spherical- harmonic 
expansion of the radiation intensity in a coordinate sys- 
tem in which the line of sight is the z direction. Note 
that we take 022 (r) to be a function of distance only, 
and not direction, consistent with our assumption that 
the quadrupole is coherent over a large patch of the sky. 
Since we use Limber's approximation below, in which an- 
gular correlations are induced only by spatial separations 
at the same distance, the variation of 022 (r) with distance 
can be included consistently. 

With the polarization written as a projection, Eq. (|g), 
along the line of sight, the angular power spectrum fol- 
lows in the flat-sky limit from Limber's equation |22]j , 

(jEE _ (jBB 

where P$ is the power spectrum of dr/dr, proportional 
to the power spectrum of the electron density n e , and the 
integral is taken up to the redshift z ro i of reionization us- 
ing the comoving differential volume given by d 2 V/dfldz. 

We now assume that the free-electron density is dis- 
tributed like the mass in the Universe and model the mass 
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distribution following the halo approach to large-scale 
structure We thus decompose the power spectrum 
into two parts, one that describes contributions from sin- 
gle halos (1-halo term) and a part that accounts for cor- 
relations between halos (2-halo term) P,p7|jlq,|2^] : 



2/1 



where 
P%(k,z) 



p(t) plh 1 p 



,, r dn(M, z) . , , .9 
dM \. I '- \n(M,z)\ 2 , 



(8) 



dM 



P%(k,z) = P^(k,z) 



dM 



dn(M, z) 
dM 



b(M, z)-n(M, z) 



(9) 



Here, (dn/dM)(M, z) is the mass function of halos pi 



1 



is the halo bias |E5| with 



and b(M, z) — , — 

v{M,z) — S c /a(M,z) the peak-height threshold with 
rms fluctuation within a top-hat filter at the virial radius 
corresponding to mass M given by a(M, z) and thresh- 
old overdensity of spherical collapse given by S c . Use- 
ful fitting functions for these quantities are tabulated in 
Ref. @. 

We define the projected scattering optical depth as 
t{9) — <jt J n e (y,6)dy where y is the line-of-sight dis- 
tance along each halo at angular distance 9 from the 
cluster center. The two-dimensional Fourier transform 
of t{9) is, 



Tl = 2ir 



9d9 t(9)J (W) , 



(10) 



where v ir corresponds to the virial radius of the halo. 
For simplicity, we model the electron distribution within 
each halo using a /3-profile and normalize it such that the 
gas mass fraction of each halo produces the global baryon 
fraction |l7) . Finally, P hn (k) is the power spectrum of 
the linear density field. We use the formulae of Ref. |27j] 
to describe the transfer function and normalize the power 
spectrum to match as — 0.9 consistent with COBE fluc- 
tuations at large scales pg| . We comment below on the 
erg dependence of our results. 
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FIG. 1. The temperature quadrupole, Cf e , as a function 
of redshift. We show both the primordial and the kinematic 
quadrupole. The bottom kinematic-quadrupole curve is for 
g(x) = 1, as appropriate for the Rayleigh- Jeans (RJ) part of 
the frequency spectrum, and the dashed and dotted curves 
are, respectively, for frequencies 150 and 220 GHz. 



quadrupole at redshift z. At redshift z — 0, with the 
power-spectrum tilt fixed at unity, COBE finds Cf^iz = 
0) = (27.5 ± 2.4 /iK) 2 |9). At higher redshifts, the mean 
quadrupole moment evolves due to the integrated Sachs- 
Wolfe effect and possibly, if the power spectrum is not 
flat, due to any scale dependence since the quadrupole 
probes smaller distances at earlier times. We calculate 
Cf e (z) following Ref. || for our fiducial ACDM cos- 
mology and show the result in Fig. [I]. 

Finally, note that Thomson scattering from cold elec- 
trons will not change the photon frequency. Thus, 
there will be no frequency dependence of the primordial 
quadrupole if we use Stokes parameters in temperature 
units — rather than the intensity units used by Ref. JTs|] — 
as we do throughout this paper. 



B. Kinematic Quadrupole 



A. Primordial Quadrupole 

The primordial quadrupole will have a coherence 
length comparable to the horizon. We thus expect that 
the amplitude of the polarization power spectrum mea- 
sured on a O(10°) patch of sky may differ by order unity 
from that on a different O(10°) patch of sky (and like- 
wise that our calculation should break down for I < 10). 
However, when averaged over the entire sky, the am- 
plitude of the power spectrum should be given by re- 
placing the quantity |et22(z)| 2 in Eq. (|t]) by its expec- 
tation value, Cf >e (z), the variance of the temperature 



As discussed in Ref. ||, an electron gas moving with 
a transverse velocity Vt relative to the CMB rest frame 
sees an 022 quadrupole moment (in a coordinate system 
in which z is along the line of sight), 

/ 4-7T 

«22=.g(x)^/-^e 2 ^, (11) 

where <p v is the orientation angle for vt on the plane of 
the sky. To obtain this result, note that in a coordinate 
system in which the z axis is aligned with the cluster's 
motion, the quadrupole dependence of the radiation tem- 
perature is g(x)v^(fi 2 - 1/3) = 3(x)w 2 (2/3) v / 47r/5F 2 ,o, 
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FIG. 2. The solid curve shows the frequency spec- 
trum of the kinematic temperature (rather than intensity) 
quadrupole. For reference, we also show the spectral depen- 
dence of the SZ thermal effect with the dotted curve. 



typical redshift of order unity, then the power-spectrum 
amplitude might generally differ by order unity from the 
power-spectrum amplitude on some other similarly sized 
patch of sky. This would be true if all the fluctuations 
were produced at one distance. However, since there will 
be 0(100) 60-Mpc coherence patches along the line of 
sight, these will tend to average out, and the polariza- 
tion power-spectrum amplitude should be roughly the 
same from one 1° degree patch to another. It is thus ap- 
propriate to replace \a,22{z) | 2 in Eq. (0) by its expectation 
value, 

(M 2 ) = %9 2 (x)(vt) = ^9 2 (x)vLs , (12) 

where we have used (vf) — {(v 2 + v 2 ) 2 ) = (8/9)vf ms , 
since (t£) = {vl) = vlJZ, = (v y ) = 3(v 2 x f = 

^msA and (v 2 x v 2 y ) = {v 2 x )(v 2 y ) = v? ms /9. We calculate 
the linear-theory rms peculiar velocity from 



v 2 (r) = 



k 2 dk / G 



2tt 2 



i) plin(M) 



(13) 



where /i is the cosine of the angle between the radiation 
direction and the z direction. However, in a coordinate 
system in which the z axis is taken to be along the line of 
sight, (^ 2 - 1/3) = -(1/3)^5^5*2,0 + v/Vl5(y 2 , 2 + 
* 2 ,- 2 )- Thus, the coefficient of y 2 2 , the component of the 
quadrupole moment that gives rise to polarization in our 
direction, is only a fraction \/3/8 of the total quadrupole 
moment. 

Unlike the primordial quadrupole, the kinematic 
quadrupole has a frequency dependence which we denote 
by g(x), where x — hv/kT CMB . Following Refs. |]l9| , p0| , 
this frequency dependence can be calculated by expand- 
ing the spectral intensity distribution of the CMB in the 
rest frame of electrons, I v oc x 3 / (e x ^( 1+v ^ — 1), in terms 
of velocity v, with 7 = (1 — v 2 )^ 1 / 2 and /1 the cosine of the 
angle between velocity and incident photon direction. We 
then obtain the frequency dependence of the quadrupole 
term, in temperature units instead of intensity units, to 
be g(x) = (x/2) coth(or/2). In Fig. |[ we show the fre- 
quency dependence of the kinematic quadrupole, which 
was neglected in Ref. |ll| 0. For reference, we also 
plot the frequency dependence of the thermal Sunyaev- 
Zel'dovich effect, 

The bulk flows associated with large-scale structure 
have coherence scales of order ~ 60 Mpc. At first, we 
might be tempted to think that if we were to survey 
some region of sky small compared with the coherence 
length ~ 1° for the peculiar velocity when projected at a 



*Note that our Figure for the frequency dependence differs 
from Fig. 4 of Ref. |h| due to our use of temperature units 
and their use of intensity units. 



where G is the linear-theory density-perturbation growth 
factor, and the overdot represents a derivative with re- 
spect to radial distance. According the halo-clustering 
model |lq] , peculiar- velocity fields are correlated over 
large distances and the nonlinear corrections to i> 2 ms 
are small. The resulting linear-theory rms kinematic 
quadrupole is also shown in Fig. |l| assuming g(x) = 1 
as relevant for Rayleigh-Jeans (RJ) part of the frequency 
spectrum when x — > 0, and for v — 150 GHz and v = 220 
GHz. 

Before moving on to discuss our results, a few words on 
how our derivations compare with those of Hu Jig ] are 
necessary. Note that our Eq. (Q) agrees with Eq. (63) 
of Ref. jl5) after converting to the logarithmic power 
spectrum, A 2 (A:) = k 3 P(k)/2ir 2 , using Limber's approx- 
imation with k = l/d,A, and substituting |a 22 (z)| 2 = 
(47r/5)Q 2 ms for the temperature quadrupole. Similarly, 
our derivation for the polarization power spectrum due 
to the kinematic quadrupole agrees with that of Ref. [pl| , 
though this comparison requires an additional step. The 
kinematic quadrupole which replaces Q 2 ms in Eq. (63) of 
Ref. @ is given in Eq. (58) as (8/45)(l - /kin)^ ms . The 
factor /kin accounts for mode coupling associated with 
the velocity field and is given by the mode-coupling in- 
tegral in the second line of Eq. (58) of Ref. |l5|]. If we 
take the limiting case of /ki n with y\ <C 1 and y 2 — » 1, 
such that the two velocity power spectra decouple from 
each other, then /k; n = 1/6 exactly and the relevant kine- 
matic quadrupole becomes (8/54)w^ ms . This is consistent 
with Eq. (12) of th is p aper. This can be verified by sub- 
stituting our Eq . jl2j ) into our Eq. (0), which recovers 
Eq. (64) of Ref. [|l5j. The fact that Hu finds numerically 
that /kin = 1/6 to a very good approximation justifies 
our assumption that the velocity and density fields are 
effectively decoupled. 
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FIG. 3. Polarization power spectra due to the rescatter- 
ing of the primordial and kinematic quadrupoles. We break 
total power in each of these cases into the 1- (lh; dotted lines) 
and 2-halo (2h; dashed lines) terms under the halo-based ap- 
proach used here. While at large angular scales correlations 
between halos dominate, at small angular scales of order few 
arcminutes and below, contributions are dominated by the 
1-halo term. For comparison, we also show primordial polar- 
ization power spectra for E and B-modes involving dominant 
scalar (E-mode) and tensor (B-mode) contributions respec- 
tively. The tensor contribution to B-modes due to inflation- 
ary gravitational waves (IGW) assumes an energy scale for 
inflation of 10 16 GeV. The long-dashed curve is the contri- 
bution to B-modes of polarization resulting from the cosmic 
shear conversion of power in E-modes, while the dot-dashed 
line labeled "Residual Lensing" represents the noise contribu- 
tion after optimally subtracting the lensing contribution using 
higher order statistics (see text for details). 

III. RESULTS AND DISCUSSION 

We summarize our results on the polarization power 
spectra in Fig. |^. Note that the secondary polarization 
discussed here contributes equally to E- and B-modes. 
While the 1-halo term dominates at arcminute angular 
scales and below, correlations between halos are impor- 
tant and determine the total effect due to secondary po- 
larization at angular scales corresponding to a few de- 
grees. The dependence of our results on the inclusion of 
the 2-halo term is consistent with the result obtained for 
the temperature power spectra from the kinetic SZ effect 
|p0| , while it is inconsistent with the thermal SZ effect, 
where contributions are dominated by the 1-halo term 
over the whole range of angular scales. The latter behav- 
ior is explained by the fact that the thermal-SZ effect is 
highly dependent on the most massive halos, while the 
kinetic SZ effect, and the secondary polarization signals 
calculated here, are irzdependent of the gas temperature 
and thus can depend on halos with a wider mass range. 

As shown in Fig. 0, the secondary E-mode polarization 



FIG. 4. The fractional contribution to polarization power 
spectra due to scattering of the primary anisotropy temper- 
ature quadrupole as a function of redshift. Here, we plot 
din Ci/d In z, for three specific values of I (f0 2 , 10 3 , and f0 4 ). 
We show the total (solid curves) as well as the 1-halo term 
(dotted curves). Note that contributions come from a broad 
range in redshift, while, with increasing I, or decreasing angu- 
lar scale, fractional contributions in the 1-halo term increase 
to higher redshifts. 

is several orders of magnitude below the E polarization 
from the surface of last scattering. The secondary polar- 
ization is therefore unlikely to be a source of confusion 
when interpreting polarization contributions to E modes. 
The amplitude of the primary effect in B modes, due to 
gravitational waves, is highly uncertain and depends on 
the energy scale of inflation |3l| . For illustration, we show 
in Fig. y the inflationary gravitational wave (IGW) sig- 
nal assuming an energy scale for inflation of E- m fi = 10 16 
GeV; the amplitude of the power spectrum scales as Ef nS . 
At large angular scales the secondary polarization is sev- 
eral orders of magnitude below the peak of this hypo- 
thetical IGW polarization signal. If the energy scale 
of inflation is lowered considerably, say to E\ n ft < 10 15 
GeV, then we might guess that the secondary polariza- 
tion could ultimately constitute a background. 

As also shown in Fig. |3|, however, there is a contri- 
bution to the B-mode power spectrum that arises from 
conversion of the primary E modes to B modes by gravi- 
tational lensing jl2| , and this is considerably larger than 
the secondary polarization. Moreover, we also show (the 
dot-dash curve) the contribution to the irreducible B- 
mode power spectrum that remains even after the lensing 
has been optimally subtracted with higher-order correla- 
tions |32|-|35||. This residual lensing power spectrum is 
considerably larger than the polarization from reioniza- 
tion. If the power spectrum is measured at a frequency 
v ~ 220 GHz, then the polarization power spectrum from 
the kinematic effect will be boosted by a factor ~ 5 from 
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FIG. 5. The fractional contribution to polarization power 
spectra due to scattering of the primary anisotropy temper- 
ature quadrupole as a function of cluster mass (in terms of 
solar mass). Here, we plot dlnCi/dlnM, for three specific 
values of I (10 2 , 10 3 , and 10 4 ). We show the total (solid line), 
the 1-halo term (dotted), and the 2-halo term (dashed). 

the g(x) — 1 power spectrum shown in Fig. |^. More- 
over, if erg = 1 (rather than the value erg = 0.9 assumed 
in Fig. then both the secondary power spectra will 
be increased for the same reasons that the temperature 
power spectra increase by a factor ~ 3.[j] Even with the 
possible frequency and as boosts, the secondary effects 
we consider here will be unlikely to be a factor for either 
gravitational-lensing or gravitational-wave studies with 
B modes. 

In Fig. ^, we show the fractional contributions to po- 
larization power spectra associated with the scattering of 
the temperature quadrupole as a function of redshift. We 
show the total and the 1-halo term for three specific val- 
ues of I corresponding to degree scales to arcminute an- 
gular scales. Though not shown here, we find consistent 
behavior for the polarization power spectrum generated 
by the scattering of the kinematic quadrupole. As shown, 
contributions come over a broad range in redshift out to 
the assumed reionization redshift of 10 with a decrease at 
highest redshifts due to the decreasing abundance of mas- 
sive halos at higher redshifts. At arcminute scales with 
I ~ 10 4 , the signal arises from halos at z > 1. A compar- 
ison of this behavior to Fig. [j] reveals that at redshifts 
greater than 1, the primary temperature quadrupole at 
the cluster positions results from a projection of the SW 



effect only. Thus, at small angular scales, scattering 
contributions come only from the quadrupole associated 
with the SW effect and not the total that includes the 
ISW effect as well. 

In Fig. ||, we show the mass dependence of the sec- 
ondary polarization signal, again using the scattering 
of the temperature quadrupole for illustration purposes. 
We show the total, the 1-halo, and the 2-halo term at 
three different values of I. While the 1-halo term is dom- 
inated by halos at the high-mass end of the mass function, 
the 2-halo term arises from a wide range in halo mass. At 
tens of arcminute scales equal contributions come from 
halo masses in the range of 10 10 to 10 14 Mq. This is 
consistent with the equivalent result for the kinetic SZ 
effect where a wide range of masses contribute. 

Note that in addition to the auto-correlation of po- 
larization, one expects secondary temperature fluctua- 
tions due to galaxy clusters to be correlated with that of 
the secondary polarization involving the E-mode. The 
temperature-polarization cross-correlation with the B- 
modes is expected to be zero based on parity considera- 
tions. We considered all combinations between secondary 
polarization and temperature anisotropies involving ther- 
mal and kinetic SZ effects and found them to be zero 
based on simple geometric arguments. 

While our simple flat-sky derivation of the secondary 
polarization anisotropy from reionization agrees with the 
all-sky approach of Ref. |p!q| , our calculational method 
complements the one used there. We use the halo model 
to describe the non-linear power spectrum of electrons 
while in Ref. |15|, electrons of the intergalactic medium 
were assumed to trace the dark-matter-density field. Our 
numerical results agree well with those of Ref. |lf| , par- 
ticularly at large angular scales where they should both 
converge to the same linear-theory calculation. Here we 
have neglected to consider the expected smoothing of the 
electron density on small scales from reheating of the 
IGM gas. However, as Ref. jl5| shows, these effects are 
easily included and reduce the power spectrum substan- 
tially only on angular scales I > 10 4 smaller than those 
we have considered here. 

While our halo-based approach is likely to be affected 
by uncertainties related to the mass function or the dis- 
tribution of electrons within halos, we expect our calcu- 
lations to accurately reflect the polarization anisotropy 
power at small angular scales. In the case of the kine- 
matic quadrupole, it is likely that we have overestimated 
the power at scales of a few degrees or more due to our 
assumption that the velocity field is coherent at such an- 
gular scales. We expect this assumption to only affect 
the 2-halo term of correlations and to result in an over- 
estimate of Ci when I is less than ~ 1000. 



^However, the boost in the polarization power spectra will 
be smaller due to the fact that much of the polarization is 
induced by electrons in smaller halos, rather than the massive 
clusters that induce the temperature fluctuation. 
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